Exact results for strongly-correlated fermions in 2+1 dimensions 
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We derive exact results for a model of strongly-interacting spinless fermions hopping on a two- 
dimensional lattice. By exploiting supersymmetry, we find the number and type of ground states 
exactly. Exploring various lattices and limits, we show how the ground states can be frustrated, 
quantum critical, or combine frustration with a Wigner crystal. We show that on generic lattices, the 
model is in an exotic "super-frustrated" state characterized by an extensive ground-state entropy. 
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Over the last few decades thousands of papers have 
been written exploring properties of itinerant-electron 
models in two spatial dimensions. Exact results, how- 
ever, for such systems at strong coupling are few and 
far between. In this paper we find the exact number and 
type of ground states in a model of spinless fermions with 
strongly-repulsive nearest- and next-nearest-neighbor in- 
teractions. The strengths of these interactions are tuned 
to give an exact supersymmetry. The supersymmetry 
not only makes our exact computations possible, but bal- 
ances competing terms in the Hamiltonian. On most lat- 
tices, this results in an exotic "super-frustrated" state. 

Our model is most transparently defined in terms of the 
supersymmetry generator Q and its hcrmitian conjugate 
Q>, which are fermionic and obey Q 2 = (Q^) 2 = 0. These 
commute with the Hamiltonian defined by H — {Q, Q^}. 
This relation is at the heart of supersymmetric quan- 
tum mechanics; a number of important results follow 
All energy eigenvalues E satisfy E > 0, because 
(s\H\s) = {.s|<5Qt| s ) + (s|Q^Q|s) cannot be negative. Any 
state with E — is therefore a ground state; it is anni- 
hilated by both Q and QK Therefore all we need to do 
to construct a many-body model with supersymmetry is 
to find a fermionic operator Q squaring to zero. 

Our degrees of freedom are spinless fermions living 
on any lattice or graph of N sites in any dimension. 
A fermion at site i is created by the operator cj with 



•4> 



6ij. The sum £\ cj squares to zero, but using 



this for Q results in a trivial Hamiltonian. The strongly- 
interacting model we will discuss was introduced in ref. 
2. The fermions have a hard core, meaning that they are 
not only forbidden to be on the same site as required by 
Fermi statistics, but are also forbidden to be on adjacent 
sites. Their creation operator is d\ — CiV<i>, where 



n a 

j next to i 



(1) 



is zero if any site next to i is occupied. A fermionic 
operator Q squaring to zero is then Q = J^dj- This 
gives a non-trivial Hamiltonian 

H = {g,Qt }= dfo + !>«>■ 

<ij> i 



The latter term has a more conventional form on a lattice 
where every site has z nearest neighbors: 



E 



V. 



<i> 



N ~ zF - 



<i> 



(3) 



where l^<i> + 1 is the number of particles adjacent to 
i, unless there are none, in which case V < i > = 0. The 
operator F — J^. d\di counts the number of fermions. So 
in addition to the hard core, the Hamiltonian includes 
a hopping term, a constant (which we keep to ensure 
ground states have E = 0), a chemical potential z, and 
repulsive interactions between fermions two sites apart. 

We use two mathematical tools to study the E = 
ground states of The first is the Witten index W 
1] . It is similar to the partition function, but includes a 
minus sign for each fermion: 

(4) 



W = tr [(-lfe-P H ] 



W is a lower bound on the number of ground states: 
it is the difference of the number of bosonic ground 
states and the number of fermionic ground states. This 
is because all energy eigenstates with E > form bo- 
son/fcrmion doublets of the same energy E but opposite 
(-l) F . The states in a doublet contribute to W with 
opposite signs and cancel, leaving only the sum of (— 1) F 
over the ground states. 

This argument shows that W is independent of j3, so 
we can evaluate it in the (3 — > limit, where every state 
contributes with weight (— 1) F . We compute this by di- 
viding the lattice into two sublattices S\ and S2; we fix 
a configuration on S\, and sum (— 1)^ for the configura- 
tions on S2. Then we sum the results over the configu- 
rations on S\. For a periodic chain with N = 3j sites, 
we take S2 to be every third site, and the remaining sites 
Si. Then the sum over configurations on any site on S2 
vanishes unless at least one of the adjacent sites on Si is 
occupied. There are only two such configurations: 

\a) = • • ■ • □ o •□ o •□ o •□ o •□ o •□ o •□ o . . . 
I7) = • • • o □ • oD • oD • oD • oD • oD • oD • . . . (5) 

where the square represents an empty site on S^- Both 
(2) I a) and I7) have / = N/3, so W — 2(—l)f, requiring 
that are at least two ground states. 
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The second tool we use is the computation of the co- 
homology Hq of the operator Q. This tool is even more 
powerful, allowing us to obtain not just a lower bound, 
but rather the precise number of ground states, and the 
fcrmion number of each. The cohomology is the vector 
space of states which are annihilated by Q but which are 
not Q of something else (in mathematical parlance, these 
states are closed but not exact) |3j. Since Q 2 = 0, any 
state which is Q of something is annihilated by Q. Two 
states \s\) and \s2) are said to be in the same cohomology 
class if \sx) = 1^2) + Q\ss) for some state IS3). 

The non-trivial cohomology classes are in one-to-one 
correspondence with the E=0 ground states 0. To see 
this, consider an energy eigenstate \E) with eigenvalue 
E > 0. If Q\E) ^ 0, then it is not in any cohomology 
class. If Q\E)=0 but H\E) ^ 0, then \E) = Q{Qp\E)/E). 
This is in the trivial cohomology class, so only the E=Q 
ground states have non-trivial cohomology. Because they 
are annihilated by both Q and , linearly indepen- 
dent E=0 ground states must be in different cohomol- 
ogy classes. Precisely, the dimension of the vector space 
of ground states (the "number" of ground states) is the 
same as that of the cohomology. Since F commutes with 
the Hamiltonian, the cohomology class and the corre- 
sponding ground state have the same fermion number. 

To illustrate these techniques, let us first generalize 
some of the one-dimensional results of ref. to a stag- 
gered (but still supersymmetric) chain. Let Q(a) = 
Qi + aQi where a is a parameter and 



N/3 



3j-2 



4 3 



N/3 

Q 2 = S Td 



^j-i- 



(0) 



Because [Q(a)] 2 = 0, the Hamiltonian {Q(a), Q^(a)} is 
supersymmetric. It deforms by multiplying the hop- 
ping term by a for hopping on or off S 2 , and multiplying 
Vu\ by a 2 when i is on 82- For a — > 00, the E = 
ground states therefore are the states where P<3j> = 
for all j. There are only two: \a) and I7) from (J5J. For 
large but finite a, there remain two ground states; for 
example \a) — l/a\a2) + 0(l/a 2 ), where \a 2 ) is the sum 
of configurations differing from \a) by shifting one parti- 
cle one site to the right. When a = 1, we know from the 
Bethe ansatz solution that there are two ground states as 
well, and that the model is gapless Q. For there 
are also two ground states, one localized on S2, the other 
with one particle for every two sites of Si. The a <C 1 
ground states spontaneously break different parity sym- 
metries than \a) and I7) do, and a = 1 is a quantum 
critical point separating the two phases. 

We find the exact number of ground states by com- 
puting the cohomology Hq by using a spectral sequence. 
A useful theorem is the "tic-tac-toe" lemma of ref. 0- 
This says that under certain conditions, the cohomology 
Hq for Q — Qi + Q2 is the same as the cohomology 
of Qi acting on the cohomology of Q 2 - In an equation, 



Hq = Hq 1 (Hq 2 ) = H\2. As with our computation of 
W, Hyi is found by first fixing the configuration on all 
sites on the sublattice 5i , and computing the cohomology 
Hq 2 . Then one computes the cohomology of Q\, acting 
not on the full space of states, but only on the classes 
in Hq 2 . A sufficient condition for the lemma to hold is 
that all non-trivial elements of H12 have the same f 2 (the 
fermion number on 52). 

We apply this theorem to the one-dimensional chain by 
using the decomposition of Q = Qi + aQ 2 given by 
Consider a single site on 5*2. If both of the adjacent S\ 
sites are empty, Hq 2 is trivial: Q2 acting on the empty 
site does not vanish, while the filled site is Q2 acting on 
the empty site. Thus Hq 2 is non-trivial only when every 
site on S 2 is forced to be empty by being adjacent to 
an occupied site. The elements of Hq 2 are just the two 
states |a) and (7) pictured above in JSJ). Both states \a) 
and I7) belong to H i2 : they are closed because Qi\a) — 
Q1I7) = 0, and not exact because there are no elements 
of Hq 2 with fi = f — 1. By the tic-tac-toe lemma, there 
must be precisely two different cohomology classes in Hq, 
and therefore exactly two ground states with / = A/3. 
Applying the same arguments to the periodic chain with 
3/ ± 1 sites and to the open chain yields in all cases 
exactly one E = ground state, except in open chains 
with 3/ + I sites, where there are none 0. 

We emphasize that for finite a, |a) and I7) are not the 
ground states themselves. A representative of a cohomol- 
ogy class is not necessarily unique, because adding Q of 
something to it does not change the class. A ground state 
is the one element in each class which is also annihilated 
by Q t. One can use this observation in principle (and 
in practice for small numbers of sites) to construct the 
exactground states from \a) and I7) as a power series in 
1/a The presence of states \a) and I7) in the ground 
states hints that the energy is lowest when particles are 
three sites apart. The chemical potential favors the cre- 
ation of more particles, but putting them two sites apart 
causes an increase in potential energy and hopping en- 
ergy. The two effects balance at an average separation of 
roughly three sites; we call this heuristic the "3-rule" . 

Having introduced the mathematical tools necessary, 
we now turn to the study of our spinless-fermion model 
on two-dimensional lattices. We find that generically, 
there is an extensive ground state entropy: the number 
of ground states increases exponentially with the size of 
the system. This indicates that the system is frustrated; 
we will explain how in the following. 

The systematics of the one-dimensional case quickly 
extend to lattices of type A3 , which are obtained from any 
lattice (or even graph) A by putting two additional sites 
on every link. Letting Si be the original sites of A and S2 
the added sites, the only states in Hq 2 and H12 are the 
two where Si is completely full, and completely empty. 
The first gives an E = ground state with / = N\ (the 
number of sites of A), while the latter gives an E = 
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state with / = L\ (the number of links in A), with a 
possible exception when L\ — N\ — 1. When A is the 
square lattice, the two ground states on A 3 have filling 
/ = N/5 and / = 2N/5. Lattices of type A 3 are the only 
two-dimensional ones we know of where the number of 
ground states does not grow with the size of the lattice. 

Another exceptional case is the octagon-square lattice 
in the first part of fig. We take L rows and M columns 
of squares (hence N — 4LM sites). Let Si consist of 
the leftmost site on every square. Then Hq 2 is trivial 
unless all the M sites on Si in a given row either all are 
occupied, or all are empty. There are 2 L — 1 such con- 
figurations which have at least one row in Si occupied. 
Because of the hard core, all the sites of S2 adjacent to an 
occupied site on Si cannot be filled, and the remaining 
sites form independent open chains of length a multiple 
of 3. Such an open chain has just one element of Hq 21 
so each of these 2 L — 1 configurations correspond to one 
element of Hq 2 and H12. Now consider the configura- 
tion where all sites on Si are empty, so that the sites 
on S2 form M periodic chains, each of length 3L. We 
showed above that Hq 2 for each of these chains has two 
independent elements. Thus Hq 2 and H12 are of dimen- 
sion 2 L + 2 — 1. Applying the tic-tac-toe lemma to this 
case is more involved, but the conclusion is that there are 
2 L + 2 M — 1 ground states, each with N/4 fermions. 




FIG. 1: Configurations obeying the 3-rule on the octagon- 
square and nonagon-triangle lattices 



We believe that on the octagon-square lattice, the 
model exhibits a combination of Wigner-crystal order 
with frustration. There are 2 L +2 M configurations of N/4 
particles which satisfy our heuristic 3-rule. 2 L of them are 
of the form displayed in fig. one can shift all the parti- 
cles in a given row without violating the rule. This illus- 
trates how frustration arises: in each row one can shift all 
the particles without violating the 3-rule. Likewise, 2 M of 
them have particles on the top or bottom of each square. 
For mysterious reasons, the state with (k x , k y )=0 is not 
a ground state, but we believe the remaining 2 L + 2 M — 1 
ordered states dominate the actual ground states. In fur- 
ther support of this claim, we analyze the discrete sym- 
metries commuting with Q. If a given element of the 
cohomology spontaneously breaks such a symmetry, the 



corresponding ground state will break it too. The ground 
states have spontaneously-broken parity symmetries like 
the Wigner crystal states in fig.H Again like the crystal, 
all but one of the 2 L — 1 ground states first considered 
spontaneously break translation symmetry in the vertical 
direction but not the horizontal; 2 — 2 of the remaining 
ground states spontaneously break translation symme- 
try in the horizontal direction. Moreover, the number 
of ground states here can be changed by requiring that 
just one site anywhere on the lattice be occupied. Con- 
sider the octagon-square lattice with one site on Si and 
its three neighbors on £2 removed; this is equivalent to 
demanding that there be a particle on this Si site. On 
this lattice there are just 2 L_1 ground states. Only in an 
ordered system should this type of change occur. 

The A3 and octagon-square lattices are exceptional: 
on all other lattices we have studied the ground-state 
entropy is extensive. In many cases (including the tri- 
angular, hexagonal and Kagome lattices), this can be 
seen by computing the Witten index W as a function of 
the size of the lattice. Employing a row-to-row trans- 
fer matrix Tm, the index for M x L unit cells is ex- 
pressed as Wl.m = tr[(Tjnf) i ]. We found by exact di- 
agonalization that the largest eigenvalues A^ ax of the 
T M here behave as A^ ax oc A M , with |A| > 1. Clearly, 
the absolute value |A| sets a lower bound on the ground- 
state entropy per lattice site. For n sites per unit cell, 
Sgs/N > \n\W LM \/{nML)~ ln|A|/n. For the triangu- 
lar lattice, S GS /N > 0.13. 

For the nonagon-triangle lattice shown in the right half 
of fig. n the extensive ground-state entropy can be ex- 
actly computed. This lattice is formed by replacing every 
other site on a hexagonal lattice with a triangle. To find 
the ground states, take Si to be the sites on the tri- 
angles, and S2 to be the remaining sites. As with the 
chain, Hq 2 vanishes unless every site in S2 is adjacent to 
an occupied site on some triangle. The non-trivial ele- 
ments of Hq 2 therefore must have precisely one particle 
per triangle, each adjacent to a different site on 52. This 
is because a triangle can have at most one particle on 
it, and (with appropriate boundary conditions) there are 
the same number of triangles as there are sites on 52- 
A typical element of Hq 2 is shown in fig. ^ One can 
think of these as "dimer" configurations on the original 
honeycomb lattice, where the dimer stretches from the 
site replaced by the triangle to the adjacent non-triangle 
site. Each close-packed hard-core dimer configuration is 
in H12, and by the tic-tac-toe lemma, it corresponds to 
a ground state. The number of such ground states e Sas 
is therefore equal to the number of such dimer coverings 
of the honeycomb lattice, which for large N is 



Sgs 
N 



1 



tt/3 



£f01n[2cos(0)] =0.16153.. 



(7) 



The frustration here clearly arises because there are many 
ways of satisfying the 3-rule. 
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For the staggered model, Q = Qi+ aQ 2 , on the nona- 
gon-triangle lattice, the dimer states are the exact ground 
states when a — ► oo. In the singular limit a = there are 
more ground states: | \I/ 2 ) with all particles on 5 2 ; 2 N / i 
ground states with one particle on each of the tri- 

angles; and additional ground states at higher fermion 
numbers as well. For < a « 1, |f 2) and e Scs — 1 of 
the |\&^ s ^ ) remain ground states, while the others develop 
energies of order a 2 . The ground-state degeneracy can be 
lifted by including terms that break the supersymmetry. 
Consider changing the intra-triangle hopping amplitude 
to 1 — e with e > 0. At a — 0, the } have energy 
E = Ne/4, so here the Wigner crystal is the unique 
E=Q ground state. For a large and e small, the leading 
piece in the effective Hamiltonian is a "flip" of 3 dimers 
around a plaquette, as in the quantum dimer model 0. 
For generic potentials this model orders [lOj, leading to 
the possibility of a quantum critical point intermediate 
between this ordered phase and the a = one. A quan- 
tum critical point indeed occurs for the chain at a = 1 
00| , and so seems possible on general lattices for a ~ 1. 

The situation on lattices with higher coordination 
number is more complicated. There are ground states 
with more particles than the 3-rule allows: the increased 
chemical potential and possibilities for hopping compen- 
sate for an increase in potential energy. For the triangle- 
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FIG. 2: "Ordered" states for the triangle-square ladder and 
the square lattice; the red squares are sublattice 2 

square ladder in fig. [2] with N — 3n + 1 sites and open 
boundary conditions, we obtained a recursion relation for 
the ground-state generating function P n (z) — trGs(^ F ): 



P n+3 (z) = 2z 2 P n (z) + z 3 P„_i(z) 



(8) 



with P = 0, Pi = z, P 2 = 2z 2 , P 3 = z 3 . This shows the 
existence of 2 ra / 3 ground states at fermion number 2JV/9, 
and also indicates additional ground states at higher fill- 
ings, up to N/4. An "ordered" state with / = 2N/9 
violating the 3-rule is given in fig. |2 but the frustration 
is evident in that there are many such states. Using the 
recursion relation, we find that the ground-state entropy 
is set by the largest solution A max of A 4 — 2A — 1 = 0, 
giving S GS /N = (In A max )/3 = 0.1110 .... 



On a square lattice of 3L x 3M sites with periodic 
boundary conditions, the situation is similar. When S 2 
consists of the red squares in fig.|21 there are two elements 
of H12 which have S 2 empty; one of them is displayed in 
fig. |3 They have 2N/9 = 2LM particles, and also vio- 
late the 3-rule. Many more ground states with different 
fermion numbers can be found by introducing various 
types of defects in this pattern, but we have not found a 
way of counting them all ||. 

Our exact results indicate that there is a new kind of 
exotic phase for itinerant fermions on a two-dimensional 
lattice with strong interactions. This "super- frustrated" 
state exhibits an extensive ground-state entropy, and oc- 
curs because supersymmetry ensures a perfect balance 
between competing terms in the Hamiltonian. Patterns 
with charge order can be distinguished in various lim- 
its and on special lattices, but the effect of (approxi- 
mate) supersymmetry in general is that defects between 
different domains come at zero (very low) energy cost. 
For example, the charge order (stripes) found for hard- 
core fermions on the square lattice |ll| becomes super- 
frustrated as the interactions and chemical potential are 
tuned to the supersymmetric point. 
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